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ABSTRACT 


This paper describes a new Rayleigh-Ritz stress analysis methodology that 
has been developed for composite panels containing cutouts. The procedure, which 
makes use of a general assumed displacement field, accommodates circular and ellipti- 
cal cutouts in biaxially loaded rectangular composite panels. Symmetric integral 
padups around the cutout can be included in the analysis. Benchmark results are 
presented to demonstrate the accuracy of the technique, and strength predictions 
based on the average stress criterion are generated and compared with experimental 
data. Finally, the stress analysis methodology is integrated into a design procedure 
for sizing integral padups around circular cutouts, and a sample problem is solved to 
illustrate its use. 


INTRODUCTION 


Cutouts of various shapes and sizes occur at numerous locations in typical 
aircraft structures. These cutouts range in complexity from simple holes designed to 
accommodate fasteners in bolted wing splices to large, reinforced openings that 
provide systems routing and access in major fuselage bulkheads. The design of these 
details in composite structures requires an accurate stress analysis technique and a 
realistic criterion to predict structural failure. 

The analysis of cutouts in orthotropic and anisotropic materials has been 
the focus of numerous research efforts over the years. Both analytical and numerical 
methods have been developed for a variety of cutouts under different loading condi- 
tions. Many of the analytical methods are based upon the use of complex stress 
functions from the theory of elasticity. Lekhnitskii (Reference 1) developed analy- 
ses for elliptical, triangular, oval, and square cutouts in infinite anisotropic 
plates under in-plane loading. DeJong (Reference 2) extended this body of work in an 
analysis of rectangular cutouts with rounded corners in infinite orthotropic plates. 
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Recently, Prasad and Shuart (Reference 3) developed a general, closed form solution 
for the moment distribution in infinite anisotropic plates containing circular or 
elliptical holes and subjected to applied bending moments. 

Numerical procedures such as finite element, boundary element, and boundary 
collocation analysis have been widely used to provide stress analysis of finite 
plates with cutouts. Hong and Crews (Reference 4) used the finite element method to 
calculate stress concentration factors around circular holes in uniaxially loaded 
finite orthotropic plates. The SY5 computer program (Reference 5) uses the finite 
element method to provide stress analysis for finite isotropic or orthotropic plates 
with circular, elliptical, square, oval, or rectangular cutouts. Generalized in- 
plane loading conditions can be accommodated, and padups around cutouts can be mod- 
eled using the first three rows of elements around the cutout. The CREPAIR computer 
program (Reference 6) can be used to analyze elliptical, slotted, and rectangular 
cutouts in finite orthotropic plates by the boundary element method. CREPAIR is 
restricted to problems involving biaxial tension/compression or in-plane shear load- 
ing conditions. The SASCJ computer code (Reference 7) uses the boundary collocation 
method to provide stress analysis for loaded fastener holes in finite orthotropic 
plates. Recently, Klang and Owen (Reference 8) have used the boundary collocation 
method to develop a stress analysis for finite anisotropic panels containing circular 
or elliptical cutouts and subjected to in-plane shear loading conditions. 

Despite the substantial body of work devoted to stress analysis of 
orthotropic and anisotropic plates with cutouts, several deficiencies remain in the 
overall analysis capability. The analytical methods, which are often well-suited for 
repeated design calculations, are usually valid only for infinite plates. The finite 
element method requires elaborate preprocessing and postprocessing routines to pro- 
vide specialized results for cutout problems. The boundary element and boundary 
collocation approaches , which are especially useful for cutouts with irregular 
shapes, are difficult to extend to problems involving reinforced cutouts. This paper 
describes a new approach to cutout stress analysis that is being developed to over- 
come some of these limitations. 

The cutout analysis methodology documented in this paper is based on the 
Rayleigh-Ritz structural analysis technique, and it makes use of a general assumed 
displacement field that can be used to treat a wide variety of cutout problems in 
finite orthotropic plates. This approach eliminates elaborate preprocessing and 
postprocessing requirements, allows for direct calculation of stress concentration 
factors at the edge of the cutout, and easily accommodates padups and reinforcements 
around the hole. As presented here, the methodology is applicable to circular and 
elliptical cutouts in biaxially loaded orthotropic plates, but it can be extended to 
include other configurations and loading conditions. 

The following sections of the paper provide a description of the Rayleigh- 
Ritz analysis methodology and a presentation of benchmark results to demonstrate the 
accuracy of the approach. Next, the analysis methodology is correlated with experi- 
mental results and used to generate strength predictions for composite tension speci- 
mens with open circular holes. Finally, a design procedure and sample problem are 
presented to show how the analysis methodology can be used in the design of integral 
padups around circular cutouts. A complete listing of the equations that are re- 
quired to implement the cutout analysis methodology is provided in the Appendix. 
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ANALYSIS METHODOLOGY 


Consider an elliptical cutout in a biaxially loaded composite panel as shown 
in Figure 1. The panel is modeled as a thin, elastic orthotropic plate under in- 
plane loading. Due to the symmetry of the cutout geometry and loading, the analysis 
is confined to the first quadrant x,y > 0 of a rectangular coordinate system with 
origin at the center of the cutout. Polar coordinates r,0 are defined with respect 
to the rectangular coordinates x,y by the usual transformation relations r - 
(x^+y^)^, $ - tan"l(y/x). 

The principle of virtual work for the loaded panel can be expressed as 

( 5 { e } T {N}dA - 5 { u ) T { t } dS - 0 (1) 

A Js 

where , 

(u) - vector of panel displacements 

{N} - vector of panel stress resultants 

{e} - vector of strains at the panel middle surface 

{t} - vector of surface tractions applied along the panel boundary 

In Equation 1, A is the panel area and S is the curve that defines the panel bounda- 
ry. The symbol, 6, represents an arbitrary variation of a quantity consistent with 

the displacement boundary conditions, which require vanishing normal displacements 
along the axes of symmetry. 

Stress analysis of the panel can be carried out using the Rayleigh-Ritz 
method. In this procedure, an assumed displacement field containing unknown parame- 
ters is developed for the panel and substituted into the principle of virtual work, 
Equation 1, to yield equations for determination of the unknown parameters. The 
stress analysis of circular or elliptical cutouts in biaxially loaded panels is based 
upon the assumed displacement field 

N M k 

u r “ 2 £ qjfk(r)cos2j0 

j-0 k— 1 J 

( 2 ) 

N M M+k 

u e - S ^ qj f k (r)sin2j 6 

k 

where u r , u^ are polar coordinate displacement components and qj are unknown parame- 
ters. (NOTE that k is used here as a superscript.) The functions f^Cr) are given by 

f k (r) - r 2(k-M)+l fork - (3) 

where the superscripted quantity is an exponent. The displacement field of Equation 
2 satisfies the symmetry conditions , which require vanishing normal displacement and 
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shear stress along the x and y axes in Figure 1. The parameters M and N in Equation 
2 can be adjusted to give the required number of series terms necessary for accurate 
solution of the problem. Substitution of Equations 2 and 3 into the principle of 
virtual work leads to a system of linear algebraic equations of the form 

[K](q} - {T} ( 4 ) 

where 

[K] = stiffness matrix 

{ q > = vector of unknown coefficients 

{T) = load vector 

Details of the derivation of Equation 4 and the subsequent calculation of panel 
stresses and strains are given in the Appendix. 

The foregoing analysis procedure can be extended to account for the presence 
of symmetric integral padups around the cutout. Consider an integral padup around a 
circular cutout as shown in Figure 2. The reinforced region, A^, has thickness tp 
and the unreinforced region, A 3 , has thickness t. These two constant thickness 
regions are joined by a tapered region, A 2 , where the thickness varies linearly with 
the radial coordinate. The stiffness matrix [K] for the reinforced panel is calcu- 
lated by carrying out the integration over the panel area, A, in three parts for 
regions A^ , A 2 , and A 3 , respectively. In Region A 2 , a linear variation of the 
extensional stiffness is assumed. Hence, 


(2) (r-R 9 ) (1) (r-R ) (3) 

Aik = * — A -jv + — A-su 

J (Ri-r 2 ) j (R 2 -R 1 ) J 


(5) 


where, R^ , R 2 are defined in Figure 2, Ajfc are the extensional stiffness for region 
i, and j ,k = 1,2, or 6 . The general procedure given in the Appendix for calculation 
of the unreinforced panel stiffness matrix can be used for the reinforced panel as 
well . 


BENCHMARK RESULTS 


The analysis methodology discussed in the previous section was implemented 
in a FORTRAN computer program, and benchmark results were generated to test its 
accuracy relative to well-established cutout stress analysis techniques. Figure 3 
shows a comparison of results for the stress concentration factor at a circular hole 
in a uniaxially loaded rectangular plate. Rayleigh-Ritz analysis results were checked 
against results generated using the anisotropic finite width correction factors 
developed by Tan (Reference 9). Both (25/50/25) quasi- isotropic (i.e., 25% 0°, 50% 
±45°, and 25% 90° plies) and (50/39/11) orthotropic laminates were considered. For 
both types of laminates, very good agreement was obtained for a wide range of hole 
diameter/plate width ratios. The Rayleigh-Ritz results shown in Figure 3 were gener- 
ated by taking N«5 , M -6 in the assumed displacement field given in Equation 2. 
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Figure 4 shows benchmark results for the stress concentration factor at an 
elliptical hole in a uniaxially loaded rectangular plate. In this comparison, 
Rayleigh-Ritz analysis results were again compared with the anisotropic finite width 
correction factors developed by Tan (Reference 9) . The aspect ratio of the ellipse 
was b/a-2.0, and a (50/39/11) orthotropic laminate was considered. Two different 
versions of Tan's solution were used. The first version, designated as the basic 
solution, was developed from an approximate stress analysis for an anisotropic plate 
containing an elliptical cutout. The second version, designated as the modified 
solution, was developed to obtain improved agreement between analysis and test data. 
Figure 4 shows good agreement between the Rayleigh-Ritz solution and Tan's basic 
solution for cases where 2b/w < 0.5. For cases where 2b/w > 0.5, the Rayleigh-Ritz 
solution lies between the two Tan solutions. The Rayleigh-Ritz results shown in 
Figure 4 were obtained by taking N^IS.M-G in the assumed displacement field given in 
Equation 2. 

Figure 5 shows benchmark results for the stress concentration factor at 
circular cutouts with symmetric padups in uniaxially loaded rectangular plates. 
Rayleigh-Ritz analysis results were compared with results for isotropic plates given 
by Peterson (Reference 10) . The Rayleigh-Ritz results were found to agree fairly 
closely with Peterson's results for padup thicknesses of up to twice the unreinforced 
panel thickness. 


EXPERIMENTAL CORRELATION AND STRENGTH PREDICTION 


The strength prediction of composite laminates with cutouts requires that an 
appropriate stress analysis procedure be combined with a suitable criterion for 
prediction of material failure. For practical applications, one of the simplest of 
these material failure criteria is the serai -empirical average stress criterion origi- 
nally suggested by Whitney and Nuismer (Reference 11). In the average stress crite- 
rion for uniaxial tension loading, laminate failure is predicted when the average 
stress over a characteristic material distance, a Q , ahead of the cutout equals the 
unnotched strength of the laminate. The stress is averaged along the line perpendic- 
ular to the applied loading that coincides with the symmetry axis of the cutout. The 
characteristic material distance, a Q , must be determined by correlation with experi- 
mental data . 

To demonstrate the use of the Rayleigh-Ritz cutout stress analysis procedure 
in an overall strength prediction methodology, the stress analysis was correlated 
with experimental data to obtain an appropriate value of a G for use in the average 
stress criterion. For this purpose, the strength data published by Tan (Reference 
12) for 8 -ply AS4/3502 (0/90/145) s open-hole tension specimens was used. For differ- 
ent hole sizes, the Rayleigh-Ritz stress analysis procedure was used to calculate 
laminate stresses in the tension specimens, and the value of a Q required to bring the 
average stress criterion prediction into agreement with the experimental failure 
stress was determined. The a Q values were generated for circular holes with diame- 
ters of 0.1, 0.25, 0.3, 0.41, and 0.61 inch. The specimen dimensions given by Tan 
(Reference 12) for each of these cases were used in the calculations. The average a 0 
value for the entire range of hole sizes and specimen dimensions was found to be a Q = 
0.09 inch. 

To verify the predictive capability of the strength prediction procedure, 
the ratio of notched to unnotched laminate strength, o^/o Q , was calculated and com- 
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pared with experimental data published by Nuismer and Whitney (Reference 13) . The 
value a o -0 . 09 inch was used in the calculations, which were carried out for 16 -ply 
T300/5208 laminates with (0/±45/90)2 s and (0/90>4 s layups and nominal hole diameters 
of 0.1, 0.3, 0.6, and 1.0 inch. The specimen dimensions given by Nuismer and Whitney 
(Reference 13) were used in each of the eight analytical predictions. Figure 6 shows 
a comparison between the analytical predictions and experimental data as a function 
of hole diameter. For the quasi- isotropic (0/±45/90)2 s laminate, the predictions 
were 10 to 20 percent conservative for the range of hole sizes considered. The 
predictions for the orthotropic (0/90)4 s laminate were approximately 20 to 30 percent 
conservative. A plot of predicted versus measured values of o^/o 0 for the eight test 
cases is shown in Figure 7. 

Based upon the limited number of data considered, the analytical strength 
prediction procedure for open hole tension specimens seems to provide a conservative 
estimate of the specimen notched strength. Improved agreement between analysis and 
experiment can be obtained by adopting a more sophisticated form of the average 
stress criterion in which the parameter, a Q , is adjusted for different layups and 
hole sizes. At the time of this writing, experimental correlation and verification 
of this semi-empirical strength prediction procedure was continuing. 


DESIGN PROCEDURE AND SAMPLE PROBLEM 


The strength prediction methodology discussed in the previous section can be 
incorporated into a design procedure for symmetric integral padups around cutouts in 
composite plates. To illustrate this procedure, consider Figure 2, which shows a 
circular cutout of radius R surrounded by a circular padup in a biaxially loaded 
composite panel. The padup has a central region of radius R^ and thickness tp which 
tapers to the unreinforced panel thickness t at radius R 2 . The object of the design 
procedure is to determine the padup layup corresponding to a specific uniaxial panel 
design strength It is assumed that the padup has the same ply percentages as 

the parent laminate . 

The design task can be carried out by (1) determining the padup region 
thickness tp corresponding to and (2) finding a practical padup layup that close- 

ly approximates the calculated tp value. The solution to the first problem can be 
obtained by iterative application of a strength prediction procedure based on the 
average stress criterion. For this purpose, consider an unreinforced panel with unit 
applied load, N xa p p (Ib/in). The average value of N x along the x-axis over the line 
segment HD < x <^D + a Q is denoted as N x av . Due to the linear elastic behavior of 
the panel, the ratio of notched strength, a jj, to unnotched strength, o Qy is 

_N _ x.ap p # (6) 

°o ^x,av 

Next, consider a panel with a padup of the type shown in Figure 2. By similar rea- 
soning, the ratio of the notched panel strength, p, for the reinforced panel to the 
unnotched panel strength is 

cr N t 

NP * x.app p ( 7 ) 

a o N* z 

x , av 
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where, N* x>av is the average stress along the x-axis over the line segment hD < x < 
^D+a 0 . The * superscript indicates that the average stress resultant N* x av is 
computed in the reinforced region that has thickness tp . Since the value of tp 
influences the operall load transfer into the padup region from the unreinforced 
portion of the panel, N* x>av is a complicated function of tp. 

In the design problem, the desired value of tp corresponds to the case where 
a NP “ ° N* For this case, combination of Equations 6 and 7 gives 


N* (t ) o* t 
x, av p N _ p 


N. 


a N 


( 8 ) 


Equation 8 must be solved iteratively to determine the padup region thickness, t p . 

For this purpose, a design function, g(tp), equal to the left-hand side of Equation 8 
is introduced. Using two initial estimates for tp, the secant method can be used to 
find tp such that g(t p ) < 5, where 6 is a sufficiently small number. 

After the padup thickness, t p , corresponding to cr*^ bas been determined, it 
remains to find a practical padup layup with the same percentage of plies as the 
parent laminate. For the case of symmetric integral padup s , material of thickness, 
(tp-t)/2, must be added to each side of the panel. Therefore, the practical padup 
layup for each side of the panel is the least number of plies that (1) maintains ply 
percentages of the parent laminate and (2) has total thickness of at least (tp-t)/2. 

A flowchart for the padup design procedure is shown in Figure 8. To 
illustrate the use of this procedure, consider a sample problem involving the panel 
geometry shown in Figure 2. The parent laminate is a 36-ply (25/50/25) AS4/3502 
panel. The geometrical parameters are L~10 inches, W==4 inches, R=0 . 5 inch, R]^=0.75 
inch, R2=l inch, and t *=0.1872 inch. The unnotched panel strength is a c =71.7 ksi and 
the panel design strength is a^^-40 ksi. A value a o -0.09 inch is assumed for the 
average stress criterion. 

For the unreinforced panel with cutout radius R=0 . 5 inch, the notched panel 
strength 27.1 ksi. Therefore, a padup is necessary if the panel is to withstand 

the design strength a^=40 ksi. Using the iterative procedure illustrated in Figure 
8 with initial thickness estimates of t p «t and tp—l.lt, the necessary value of t p was 
found to be tp«=0.2996 inch. The practical padup corresponding to this value of t p 
has 12 plies (three 0°, six ±45°, three 90°) of material on each side of the panel. 

For the practical padup, t p ~0.312 inch. A strength analysis for this value 
of tp gives a reinforced notched panel strength of ajgp-41.4 ksi. Since this value of 
^NP g^ ves a positive margin of safety with respect to the design strength a^jj-40 ksi, 
the padup design procedure has achieved its purpose. 


SUMMARY AND FUTURE PLANS 


In the previous sections, a Rayleigh-Ritz stress analysis methodology for 
composite panels with cutouts has been described, and its applications to panel 
strength prediction and design of symmetric integral padups around cutouts have been 
illustrated. In its current state of development, this methodology can be applied to 
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circular or elliptical cutouts in rectangular panels under biaxial tension/ compres- 
sion loading. The advantages of the procedure are its ability to provide accurate 
numerical stress analysis solutions with a minimum of input data, and its capability 
to accommodate cutout reinforcements, such as padups . Its disadvantage is its re- 
stricted applicability, which currently includes only a limited number of cutout 
geometries and panel loading conditions. 

Future development efforts for the Rayleigh-Ritz cutout stress analysis 
methodology will be focused on the incorporation of substructural effects into the 
analysis, and the extension of the procedure to more complex cutout geometries and 
loading conditions. The effect of integral stiffeners can be included in the analy- 
sis by assuming strain compatibility between the stiffeners and the panel. The 
strain in the stiffeners can then be calculated using the assumed displacement field 
given in Equations 2 and 3, and virtual strain energy terms for the stiffener ele- 
ments can be added to the left-hand side of Equation 1. By this approach, the effect 
of any number of stiffener elements can be included in the analysis without increas- 
ing the number of degrees of freedom in the solution process. An effort to implement 
the substructural effects analysis for the common case of a "picture - frame” stiffen- 
ing arrangement around a cutout is currently underway. 

Future analysis efforts will also include attempts to extend the existing 
methodology to include more complex cutout shapes and loading conditions. The analy- 
sis of complex cutout shapes, such as slotted holes and rectangles with rounded 
corners, will require some modification of the basic analysis procedure presented 
here. One promising approach that will be investigated is use of the penalty func- 
tion method (Reference 14) in conjunction with the assumed displacement field of 
Equations 2 and 3. This approach, which involves constructing penalty functions from 
specified constraint conditions, will make it possible to closely satisfy stress -free 
boundary conditions along the cutout perimeter without introducing additional degrees 
of freedom. Further modifications will be required to enable the stress analysis 
procedure to accommodate more complex loading conditions. For in-plane shear load- 
ing, a different assumed displacement field will be constructed to exploit symmetries 
in the deformation of the panel under this type of loading. This displacement field 
will be implemented into the analysis using the procedures illustrated in the Appen- 
dix. Stress analysis results for problems involving generalized in-plane loading can 
then be obtained by superposing results obtained from the separate analyses for 
biaxial and shear loading. 

In future work, a concurrent approach to the development of design 
procedures and enhanced stress analysis methodologies will be pursued. The initial 
portion of this work, which consists of a design procedure for symmetric integral 
padups around circular holes in uniaxially loaded panels, was illustrated in this 
paper. A generalization of this design procedure to permit sizing of cutout rein- 
forcements and stiffening elements will be developed as the extensions to the stress 
analysis methodology are completed. 
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APPENDIX 


For problem solution, it is convenient to represent the assumed panel 
displacements given in Equation 2 in the matrix form 


(u) 


u r 


j!o [u J llq J' 


(A-l) 


where (qj) is a vector of dimension 2M containing unknown parameters and [Uj ] is a 
2 x 2M matrix of the form 


[Uj] 


f^cos2jS . . . %cos2j^ 0 


. 0 


0 


0 fisi^j* . . . f M sin2jtf 


(A- 2) 


The strain vector { e } is obtained by applying the strain-displacement, relations for 
polar coordinates : 


r ■> 

€ r 


' S_ 

dr 

0 

€ o 

► K 

1 

1 2- 

U 


r 

r 36 

, Trfi . 


1 d- 

~ <3/5 

L. - 1 

Q w J 


(A-3) 


Substitution of 

Equations A-l and A- 2 

into 

Equation A-3 

leads 

to a relation of the 

form 

N 






U) - 

jEo 1 H j |lq J' 




(A-4) 

where [Hj ] is a 

3 x 

2M matrix whose entries 

are given by 





f' cos2j 6 

for 

m - 1 , . . . 

,M 


H j ,1m 

«* < 

i m 

i 




(A-5) 


[ 0 

for 

m - M+l , . . . 

2M 




r lf n cos2jtf 

for 

m - 1 , . . . 

,M 


H j ,2m 


1 r 




(A-6) 


1 2jif m . M c°s2jtf 
r 

for 

m - M+l , . . . 

2M 




-2j — f m sin2j 6 

for 

m - 1 , . . . 

,M 


H j ,3m 

* 4 

r 

L (f’ - lf m )sin2j* 

for 

m - M+l , . . . 

2M 

(A-7) 


m r 


The superposed prime denotes differentiation with respect to r. The stress resultant 
vector (N) is obtained by application of the stress-strain relations. Let E x , Ey, 
G X y, ^xy the orthotropic elastic constants for the panel. Let 
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l 



An - 


E t 
x 


i-^xy^yx 


a 12 


v E t 

xy y 

I'^xy^yx 


(A- 8 ) 


a 22 


E t 
-Y-. ■ 


1 - I'xy^yx 


A 66 “ ^ X yt 


where u yx “ ^xy(Ey/E x ) and t is the panel thickness. The stress resultant vector is 
given by 


where 



' N r 


(N) - - 

N* 



. N rtf . 



An A 12 a 16 
a 12 a 22 a 26 
a 16 a 26 a 66 


l IrB J 


(A- 9 ) 


A 11 “ A ll cos ^ + 2(A]^2 + 2A^6) sin^0 cos^ + A22 s ^ n ^^ 

a 12 “ ( A i i + A 22 - 4 A 55 ) sin^9cos^9 + A 12( cos ^0 + sin^0) 

A^^ - - ( A n - A 12 - 2 A 55 )sin^cos^^ - (A ^2 “ A 22 + 2 A 5 g) sin^flcos# 

a 22 ” A^sin^S + 2(A^2 + 2Agg) sin^flcos^fl + A22 COS ^0 

a 26 “ “( A ll “ a 12 ' 2Agg)sin^0cos0 - ( A 12 “ a 22 + 2Agg) sinflcos^fl 

A 55 =* (A^i + A 22 “ 2A]^2 “ 2 A^ 5 ) sin^^cos^^ + Ag^sin^# + cos^tf) 

Equation A- 9 can be written in the compact matrix form 

(N) - [A] { e } 


(A- 10) 


(A- 11) 


where [A] is the 3x3 matrix in Equation A- 9. 


For the rectangular panel of length L and width W shown in Figure 1, the r 
coordinate defining the boundary is a function of 9 given by 


Rb(*) 


L 

2 cos 0 


for 0 < 9 < <f> 


W 

2 sin 0 


for 4> < 8 < - 
2 


(A- 12) 


where <f> - tan‘^(W/L) . The r coordinate defining the edge of the elliptical cutout is 
also a function of 9 given by 

R c (0) - ^ (A- 13) 
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where a,b are the dimensions of the ellipse semiaxes. The surface traction vector 
along the panel boundary is 


{t} 


t r 


c , e 


< ► 

1 1 N x (l+cos20 ) 

1 N v ( 1 - cos20 ) 

2 y 

f 'I N x sin20 
, 2 

I N v sin20 

2 y 


for 0 < 0 < <f> 

for <f> < 8 < - 

2 

for 0 < 8 < <f> 

for <f> < 8 < — 
" “2 


(A- 14) 


(A- 15) 


(A- 16) 


Introducing Equations A-l, A-4, A- 11, and A- 14 into Equation 1 leads to a matrix 
equation of the form 

[K] { q } - (T) , (A- 17 ) 


where [K] is an M(2N+1) x M(2N+1) stiffness matrix, {q} is an M(2N+1) x 1 vector of 
unknown displacement field parameters and { T } is an M(2N+1) x 1 load vector. 

The stiffness matrix [K] can be represented in the partitioned form 


[Koo] [Koi] ... . . [Ron] 

[ K iol t K ii 1 • • • ■ • [ k in1 


[K] - 


(A- 18) 


L I k NQ] [%l] • • • • • [ k nn1 J 


where the partition [KqoI is an M x M matrix, the partitions [Kqi] , . . . , [Kqn] are M x 
2M matrices, the partitions [Kiq] , • . . , [% ol are 2M x M matrices, and the remaining 
partitions are 2M x 2M matrices. The elements of a general partition [ Kj y- ] are given 
by the relations 


K j k , mn 


*7T 

2 UAn/9 mn + A 12 (6 mn + 6 n m ) + A 22 ^ mn ] cos2j 0cos2k0 

. 0 


- 2k(A 16 5 nra + A 26 ^ mn )cos2j0sin2k0 

- 2j(A 16 6 mn + A 2 6^ mn )sin2j 0cos2k0 
+ 4jkA66^ mn sin2j Ssin2kS }dS 

for 1 < m < M, l<n<M 


(A-19) 
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where , 


^ j k , ran 


^ j k , mn 


^ j k , mn 


*7T 

2 t 2k ( A 12 5 n-M,m + A22^m,n-M) cos2 j0 cos2R 0 

Jo 


+ ( A 16 A m,n-M + A 26^m , n-M) cos2 j 0sin2k0 

- 4jkA26^ m) n-M sin2 j^cos2k^ 

- 2 J A 66 , ?m > n-M sin2 j^ si -n 2 k^ ]A9 

for 1 < m < M, M+l < n < 2M 


'n 

2 t 2 j ( A 12 5 m-M,n + A 22^m-M,n) cos2 j^ cos2 k^ 
. 0 


- 4jkA26<£ m _M n cos2j0sin2k0 

+ ( A 16^n,m-M + A 26^n,m-M) sin2 j 0cos2k0 
^kA^ fat) j-j ^ in2 j 8 s in2k0 ] d.0 

for M+l < m < 2M , 1 < n < M 

'n 

2 [4jk A 22<^ m -M ,n-M cos2 j 0cos2k0 

. 0 

+ 2 J A 26 , 7m-M ) n-M cos2 j0 sin2 k0 
+ 2kA26r/ n . M>m . M sin2j$cos2kS 
+ A 66 w m-M,n-M sin2 j 0sin2k0 ]d0 

for M+l < m < 2M , M+l < n < 2M 


f R b , , 

£mn ” f m f n rdr 

J R c 

f R b 

5 mn “ f m f n dr 

J R c 

f R b 

^mn “ if m f n dr 

J R c r 

f R b , , 

-^mn “ ^m(fn — fn) rdr ” $mn " ^nm 

JR^ r 


(A-20) 

(A-21) 

(A-22) 

(A- 23) 
(A-24) 
(A" 25) 
(A-26) 
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^mn ” 


(A-27) 


w mn 


R b « 

An^n ' ” ^mn * Ann 

J R c r 

>R b 

(fm ' — Ann ‘ ^mn ' ^nm + ^mn 

p r 


(A- 28) 


for 1 < m,n < M. In Equations A-23 through A-31, the prime denotes differentiation 
with respect to r. Using the functions defined in Equation 3, the integrals in 
Equations A-23 through A- 25 can be evaluated to obtain 


^mn 


r2(m-MUn f2(n-M)+ll 
[ 2 (m+n) - 4M+2] 


[2 (m+n) -4M+2] 
• R b 

L 


[2 (m+n) -4M+2] 1 
R c \ 


for m+n * 2M-1 , 


(A-29) 


[2(m-M)+l] [2(n-M)+l] ln(R b /R c ) 

for m+n — 2M-1 


S 


mn 


^mn 


[9/n MU11 f [ 2 (m+n) -4M+2 ] 
I 2(n-M)+l I J R , 

[2 (m+n) - 4M+2 ] [ 


[2 (m+n) -4M+2] 
R c 


for m+n 2M-1 , 


[ 2(n-M)+l ] ln(R b /R c ) 


for m+n — 2M-1 


1 

[2 (m+n) - 4M+2] 


f [2 (m+n) -4M+2] 
\ R b 


[ 2 (m+n) -4M+2 ] 1 
R c \ 


for m+n * 2M-1 , 


(A- 30) 


(A-31) 


ln(R b /R c ) 


for m+n — 2M-1 


The quantities R b and R c are functions of 9 as defined in Equations A-12 and A-13. 
The integrals in Equations A- 19 through A- 22 can be evaluated numerically using 
Simpson's rule. 

The load vector (T) in Equation A- 17 can be represented in the partitioned 

form 
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(T) 


' {T 0 J ' 
{T X } 

(Tj) 

. (T N ) j 


(A-32) 


where (Tg) is an M x 1 vector and the remaining partitions are 2M x 1 vectors. The 
elements of a general partition (Tj } are given by the relations 


L J ,m 


_ 1 


± N V L 


f m (R b )cos2j*d0 


+ -i N„W 


^ f m (R b ) cos 2j for m - 1 M 

4> 




- - n x l 
2 




f m -M(^b) s i n 2j ^tan^dS 


+ i N V W 


^ fm-M^b) s ^- n 2j ^cot^d^ for m — M+l, . 


. ,2M 


(A- 33) 


(A- 34) 


These integrals can be evaluated numerically using Simpson's rule. Finally, the 
vector of unknown displacement field parameters { q ) in Equation A- 17 can be written 
in the partitioned form 


<q> 


{qo> ' 

(qi) 


Un* J 


(A- 35) 


where ( qo } is an M x 1 vector and the remaining partitions are 2M xl vectors. 


Solution of the system of linear algebraic equations represented by Equation 
A-17 yields the vector of displacement field parameters (q). Displacements, strains, 
and stress resultants at any point in the cutout panel can be calculated by evaluat- 
ing Equations A-l, A-4, and A- 9, respectively. 
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Figure 3. Benchmark results for finite structure circular cutout analysis. 



Figure 4. Benchmark results for finite structure elliptical cutout analysis. 
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Figure 5. Benchmark results for uniform thickness circular padups. 
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Analytical Predictions 
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Figure 6. Analytical strength predictions and experimental data for different hole sizes 
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Figure 7. Analytical strength predictions versus experimental data. 





Figure 8. Flowchart for Padup Design Procedure. 
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